He-LiF surface potential from fast atom diffraction under grazing incidence. 
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Diffraction patterns produced by grazing scattering of fast atoms from insulator surfaces are 
used to examine the atom-surface interaction. The method is applied to He atoms colliding with a 
LiF(OOl) surface along axial crystallographic channels. The projectile-surface potential is obtained 
from an accurate DFT calculation, which includes polarization and surface relaxation. For the 
description of the collision process we employ the surface eikonal approximation, which takes into 
account quantum interference between different projectile paths. The dependence of projectile spec- 
tra on the parallel and perpendicular incident energies is experimentally and theoretically analyzed, 
determining the range of applicability of the proposed model. 

PACS numbers: 79.20.Rf, 79.60.Bm, 34.20.Cf. 



I. INTRODUCTION 

Diffraction of thermal atoms from crystal surfaces has 
been extensivel y s tudied over the years [l], H, H, 0, H, 
H, 0, B H EH HH , becoming a common tool for sur- 
face analysis. Recently new experiments [HI, [l3|, Q, HH 
have shown interference effects also for grazing scat- 
tering of fast atoms from surfaces, where classical me- 
chanics was supposed to be adequate. This unexpected 
diffraction phenomenon was found to be very sensitive 
to the description of the projectile-surface interaction 

M M TH El which opens the way for a method 

to probe surface potentials with high accuracy. 

The aim of this work is to find out to what extent sur- 
face potentials derived from statc-of-art ab-initio meth- 
ods are capable of reproducing experimental diffraction 
patterns for grazing scattering of swift He atoms from 
a LiF(OOl) surface. The He-LiF surface interaction is 
here derived by using the Siesta pol ] implementation 
of the density- functional theory (DFT), which is a self- 
consistent method for performing first-principles calcu- 
lations on systems with a large number of atoms. This 
DFT method has been successfully used to study a va- 
riety of nanoscale problems [21j . In order to describe 
the interference process, we employ a distorted-wave 
model - the surface eikonal approximation [l6| - using 
the eikonal wave function to represent the clastic collision 
with the surface, while the motion of the fast projectile 
is described classically by considering axially channeled 
trajectories for different initial conditions. The surface 
eikonal approximation is valid for small de Broglie wave- 
lengths of incident atoms, as considered here, which are 
several orders of magnitude smaller than the interatomic 
distances in the crystal. This method was shown to pro- 
vide an adequate description of the interference effects 
for atoms colliding with insulator surfaces under axial 



surface channeling [22J. 

Eikonal projectile angular distributions derived from 
using the DFT surface potential arc compared with the 
experiment for different energies of incident projectiles. 
From this comparison we deduce the validity range of 
the potential model, which involves polarization and sur- 
face rumpling. The paper is organized as follows. The 
experimental method and the theoretical formalism are 
summarized in Sec. II and III, respectively. Results are 
presented and discussed in Sec. IV, and in Sec. V we 
outline our conclusions. Atomic units (a.u.) are used 
unless otherwise stated. 



II. EXPERIMENTAL METHOD 

In our experiments we have scattered neutral 3 He and 
4 He atoms with kinetic energies Ei ranging from 0.3 
kcV to 25 keV from a clean and flat LiF(OOl) surface 
at room temperature under grazing angles of incidence 
0.4 < 9i < 1.5 deg. Fast Hc + ion beams were produced in 
a 10 GHz electron cyclotron resonance (ECR) ion source 
(Nanogan Pantechnique, Caen, France). The neutraliza- 
tion of the He + ions was achieved via charge transfer in a 
gas cell mounted in the beam line of the accelerator oper- 
ating with He gas and subsequent deflection of remaining 
charged fraction by an electric field. A base pressure of 
some 10 -11 mbar was achieved in our UHV chamber by a 
turbomolccular pump in series with a titanium sublima- 
tion pump, where the pressure gradient with respect to 
the beam line of the accelerator was maintained by two 
differential pumping stages. Pairs of slits at both ends of 
these stages were used for the collimation of the incident 
beam to a divergence of < 0.03°. This high collimation is 
necessary for diffraction in order to maintain the degree 
of coherence in the scattering process from LiF(OOl) . 
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The LiF(OOl) surface was prepared by cycles of graz- 
ing sputtering with 25 keV Ar + ions at 250° C where the 
ionic conductivity of LiF is sufficiently enhanced in order 
to avoided macroscopic charging up and subsequent an- 
nealing to temperatures of about 350°C. The scattering 
experiments were performed in the regime of axial sur- 
face channeling, i.c the azimuthal setting of the surface 
plane was chosen so that the direction of the incident 
beam was parallel with atomic strings along low indexed 
directions in the surface plane. 

2D angular distributions of scattered projectiles were 
recorded by means of a commercially available position- 
sensitive multi-channelplate detector (MCP) with a 
delay-line anode (DLD40, Roentdek Handels GmbH) lo- 
cated 66 cm behind the target. This provides a sim- 
ple and very efficient procedure for recording data where 
complete diffraction patterns as shown below can be 
recorded in a time of about minutes. Since only about 
10 4 He atoms per second hit the target surface, fast atom 
diffraction is non-destructive and can be applied in stud- 
ies on insulator surfaces (neutral projectiles) [lH. Li. as 
well as adsorption phenomena at metal surfaces 15l . 18| . 

Since aside from the absolute angular positions of 
diffraction spots, their relative intensities are important 
here, one has to carefully correct the recorded diffraction 
patterns with respect to inhomogencities in the detection 
efficiency across the sensitive area of the MCP. This cor- 
rection is performed by means of a wobbling technique 
where the projectile beam is scanned across the MCP 
active area with two orthogonal oriented electric fields 
using frequencies in the kHz domain. 



III. THEORETICAL MODEL 

When an atomic projectile (P) impinges under grazing 
incidence on a crystal surface (S) the T- matrix element 
associated with the elastic scattering can be defined in 
terms of the scattering state of the projectile, , as 



where 



l if 
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v(Kp(t))= / dt' Vsp{Hp{t'))+4> M (3) 



is the eikonal-Maslov phase, which depends on the clas- 
sical position of the incident atom at a given time 
t. This phase includes the Maslov correction term <fi M = 
un/2 that takes into account the phase change suffered 
by the projectile at tur ning points, with v the Maslov 
index defined as in Ref . [261] . By inserting Eq. in 
Eq. OJ), a fter some algebra the cikonal transition matrix 
reads [16| 



T 



(eik) 
if 



dR os a,if(R os ), 



(4) 



where R os determines the initial position of the projectile 
on the surface plane and 
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is the transition amplitude associated with the classi- 
cal path Kp(R os ,t). In Eq. © Q = Kf - K t is the 
projectile momentum transfer and v z (1Zp) denotes the 
component of the projectile velocity perpendicular to the 
surface plane, with z directly along the surface normal, 
towards the vacuum region. 

The differential probability, per unit of surface area, 
for elastic scattering with final momentum Kf in the di- 
rection of the solid angle f2/ = (6*/ , V 9 / ) l& obtained from 



rfi(eik) 
1 if 



where T, 
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de- 



Eq. dU as dP/dflf = (2tt) 4 5 

notes the eikonal T-matrix element, normalized per unit 
area, rap is the projectile mass, and Of and ip^ are the 
final polar and azimuthal angles, respectively, with ipj 
measured with respect to the x axis, along the incidence 
direction in the surface plane. Details are given in Rcfs. 

MM- 



where Vsp is the surface-projectile interaction, 
Rp denotes the position of the center of mass of the 
incident atom, and $j(Rp) = (2tt)~ 3 ^ 2 cxp(iKj ■ Rp), 
with j = i(f), being the initial (final) unperturbed 
wave function and Kuf\ the initial (final) projectile 
momentum. Taking into account that in the range of 
impact energies the dc Broglie wavelength of the incident 
projectile, A = 2n/Ki, is sufficiently short compared to 
the characteristic distance of the surface potential, we 
approximate the scattering state ^>f by means of the 
eikonal-Maslov wave function [25| as follows 

*+(£p) ^ xT k)+ {K P ) - $ l (7e P )cxp(- l 7 7 (7?p)), (2) 



A. Projectile-surface interaction 

The surface potential was determined by performing 
first-principles calculations for the LiF(OOl) surface. We 
used the Siesta [2p| implementation of density-functional 
theory (DFT) within the Local Density Approximation 
(LDA) to obtain the effective interaction potential 
between a He atom and a slab of 10 atomic planes of 
LiF. Periodic boundary conditions were used in the (001) 
plane, and in order to prevent the interaction between 
images of the He atoms a y/{2) x y/{2) supercell was 
considered in the (001) plane, giving a total of 112+1 
atoms in the simulation box. 
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We model the core electrons with norm-conserving 
pseudopotentials of the Troullier-Martins type (24|, de- 
scribing valence electrons with numerical double-zeta po- 
larized atomic orbitals as the basis set. We explicitly in- 
cluded semicore states for Li (Is 2 ) in the valence band, 
and added two extra layers of "ghost orbitals" at the sur- 
face to increase the basis set description of the surface 
electronic wave-functions. 

The structure was optimized until the forces on all 
atoms were smaller than 0.03 eV/A. The slight underes- 
timation for the in-plane lattice constant obtained in our 
slab geometry (3.94A, experimental 4.02A Q) is typical 
of LDA calculations. Different relaxation of the Li and 
F atoms at the surface called rumpling is apparent in 
the first two atomic layers, with F atoms in the surface 
plane being slightly pushed out and Li atoms slightly de- 
pressed by a distance d\ measured with respect to the un- 
reconstructed surface. For the topmost atomic layer we 
obtained d\ = 0.046 a.u., while the displacement corre- 
sponding to the second layer is substantially smaller and 
in opposite direction (d,2 = —0.010 a.u.). These values 
are slightly higher than those reported from a IV-LEED 
analysis [29j . 

The surface potential V(z) for a given position (x,y) 
on the surface plane was obtained from calculations of 
the total energy of the system composed by the LiF slab 
and the He atom placed a distance z from the last atomic 
layer of the slab (taken as the average between Li and F 
positions). The standard correction [3(| due to BSSE 
(Basis Set Superposition Errors) was considered for the 
computed V(z): 

V(z) = £[LiF+He(z)] - £[Hc(z)] LiF - £[LiF] Hc(z \ (6) 

where i?[A] x denotes the energy of the system A consid- 
ering not only the basis orbitals of A, but also those that 
correspond to the subsystem X. Once V(z) is known for 
selected 9 in-plane high-symmetry positions in a mesh 
of z points, an interpolation scheme is used to derive 
Vsp(R) at any point in the vacuum region. 

IV. RESULTS 

Experimental angular distributions of He projectiles 
elastically scattered from a LiF(OOl) surface under ax- 
ial surface channeling conditions are presented here as a 
benchmark for the DFT projectile-surface potential. 

First, we analyze the dependence of final projectile 
distributions on the incidence energy Ei, splitted into 
two terms £7, = Em + Eij_, where Eq = Ei cos 2 8 
(Eij_ = E sin 2 9i) is associated with the component of 
the initial velocity parallel (perpendicular) to the axial 
channel, with 9i the incidence angle measured with re- 
spect to the surface plane. In Fig. [I] we show diffraction 
patterns for 4 He atoms impinging along the (110) direc- 
tion with two different energies - Ei = 2.2 keV and 7.5 
kcV - but with the same perpendicular energy, Ei± =1.04 



FIG. 1: Two dimensional intensity distributions, as recorded 
with a position sensitive detector, for 4 He atoms scattered 
from LiF(OOl) along a (110) direction with two different pro- 
jectile energies, Ei~2.2 keV and 7.5 keV, but with the same 
perpendicular energy (Ei± — 1.04 eV). Color code: Red=high, 
blue=low intensity. Positions of rainbow angles Q r b are in- 
dicated by straight lines. Black dots represent theoretical 
positions of maxima from the surface eikonal model. 



cV. In both cases, the distribution for scattered projec- 
tile lies inside an annulus of radius 6i , presenting maxima 
symmetrically placed with respect to the incidence direc- 
tion (i.e. Lpj — 0). The high intensity for the outermost 
peaks (marked with lines) is due to rainbow scattering 
under rainbow angle r b at maximal deflection, which 
can be explained classically. Between the outer rainbow 
peaks further peaks show up, which can be explained as 
quantum mechanical diffraction effects in close analogy 
to the origin of supernumerary rainbows [l4j . Super- 
numerary rainbows originate from quantum interference 
between projectiles that follow different classical path- 
ways with the same final momentum. The order m of a 
supernumerary corresponds to the multiple of A in path 
length difference for constructive interference. The black 
dots in Fig. [1] represent the theoretical positions of the 
maxima from the eikonal model, which closely agree with 
experimental data. 

In Fig. [5] we show the intensity inside the annulus of 
radius Qi from Fig. [T]as a function of the deflection angle 
O, defined as O = arctan((^ /Of). Position and num- 
ber of the supernumerary maxima are independent of Ei 
at the same Ei±. Similar structures are predicted by 
the eikonal model, although in the vicinity of the clas- 
sical rainbow angle Q r b the relative intensity is overes- 
timated. This is because the eikonal model is a semi- 
classical method based on classically calculated trajecto- 
ries, showing a sharp maximum at the classical rainbow, 
where intensity increases sharply for O — > r j, and is 
zero for > Q r b- In a more elaborate quantum treat- 
ment the classical rainbow peak will be replaced by a 
smoother maximum for m = at < Q r b, with de- 
caying intensity on the dark side of the classical rainbow 
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FIG. 2: Projected intensities inside annuluses in Fig. [T] as 
a function of the deflection angle O, for Ei=7.5 keV (gray 
circles) and 2.2 keV (blue squares) and corresponding differ- 
ential probabilities derived from the surface eikonal approach 
(black dashed and red full curves), m denotes the order of 
the supernumerary rainbow. 



account the complete corrugation on the surface plane, 
without averaging the projectile-surface potential along 
the incidence direction, the eikonal projectile distribution 
is practically unaffected by the modulation of the poten- 
tial along the channel. Thus the differential probability 
dP/dQ is independent of E^i for a given perpendicular 
energy. On the other hand, the positions of supernu- 
merary maxima are found to be extremely sensitive to 
the shape of the surface potential across the channel, es- 
pecially for higher orders m, which correspond to small 
deflection angles [U, EE HI • From Fig. [3] we show 
that the DFT surface potential reproduces fairly well the 
angular positions of the supernumeraries. 

Different perpendicular energies Ei± probe a different 
z-range of Vsp. To investigate in detail the atom-surface 
potential across the (110) channel, in Fig. 3] we plot the 
angular positions of maxima of the experimental distri- 
bution as a function of Ei± ranging from 0.03 to 3 eV. 
In this case 3 He isotopes are used as projectiles. For low 
perpendicular energies, i.e. Ei± < 1 eV, the experimental 
spectra show maxima at Bragg angles 0„, which fulfill 
the condition 



d sin 0„ = n\_ 



(7) 



FIG. 3: Deflection angles O corresponding to maxima of an- 
gular distributions, as a function of the projectile energy Ei, 
for He atoms scattered from LiF(OOl) along the direction 
(110). The perpendicular energy is kept as a constant (Ei± 
=1.04 eV). Circles, experimental data; curves, quantum rain- 
bow m = and supernmerary rainbows m = 1 to 4 derived 
within the surface eikonal approximation. 



G > Q r b [3ll . l32j . The maximum in = is the "quantum 
surface rainbow" [33[. 

Just like in the experiment, eikonal patterns as a func- 
tion of the deflection angle are independent of Ei at the 
same Ei±. This becomes more evident in Fig. [3] where 
deflection angles corresponding to supernumerary rain- 
bows are plotted as a function of the total energy Ei while 
keeping Ei± constant. Even though Vsp(R) takes into 



d being the width of the channel, n the diffraction or- 
der, and Aj_ = 2ir/Ki Z the de Broglie wavelength asso- 
ciated with the perpendicular motion. As discussed in 
Ref. 0, [3 , interference patterns for grazing scattering 
stem from two different mechanisms. The first one, as- 
sociated with the supernumerary rainbows, is produced 
by the interference of trajectories whose initial positions 
R os differ by a distance smaller than d, carrying informa- 
tion on the shape of the interaction potential across the 
channel. The second one originates from the interference 
of trajectories whose initial positions R os are separated 
by the spacial lattice periodicity d resulting in "Bragg 
peaks" providing information on the spacing between sur- 
face atoms. Whether a Bragg peak shows intensity or 
not depends on the position of the supernumeraries. The 
Bragg peak of order n which is closest to the angular 
position of a supernumerary of order m is intense. Since 
the Bragg angles 0„ decrease with Ei± while the angular 
position of the supernumeraries increases, the order n of 
the intense Bragg peak increases successively. 

Signatures of both interference processes can be ob- 
served in the simulated spectrum also. In Fig. [5] the 
eikonal probability dP/dQ is plotted as a function of the 
deflection angle for E.- L ±_ = 0.5 eV. From Eq. ^ when 
the initial projectile position R os is integrated over a unit 
cell, supernumerary maxima are only present in the an- 
gular projectile distribution. But when the integration 
area is extended to include the first order nearest neigh- 
bor target ions, the eikonal spectrum displays internal 
structures in the supernumerary maxima, which are due 
to resolved Bragg peaks. 

Since we are interested in studying the modulation of 
the potential inside an unit cell, we have plotted in Fig. 
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FIG. 4: Similar to Fig. [3] but as function of the perpendicular 
energy Ei±, for 3 He atoms scattered from LiF(OOl) along the 
direction (110). Full colored curves, quantum rainbow m — 
and supernmerary rainbows m = 1 to 4 derived within the 
surface eikonal approximation based on the present DFT po- 
tential; dashed curves, positions of supernumerary rainbows 
in the hard wall approximation (Eq. (jSJ)) using the effective 
corrugation of the DFT potential across the (110) channel; 
grey curves, theoretical positions of maxima from Bragg con- 
dition (Eq. ©). 



[4] eikonal curves corresponding to the center of supernu- 
merary maxima, neglecting the Bragg interference that 
appears as a superimposed structure at low perpendic- 
ular energies. The eikonal curves obtained by using the 
DFT atom-surface interaction follow closely the experi- 
mental results for supernumeraries m = 1 to m = 4 with 
only slight deviations for the quantum rainbow m = at 
low Ea_. Supernumerary rainbows are sensitive to the 
corrugation of the equipotcntial surfaces but not to their 
positions. Each potential with the same effective corru- 
gation along the (110) direction would agree as well. In 
order to show that the potential is unique, one has to 
compare the theoretical results with the experiment for 
a further (different) channeling direction. The potential 
would be unique only if the same potential describes the 
supernumeraries here as well. 

In Fig. [6] we show a comparison of experimental su- 
pernumerary rainbows with eikonal results for scattering 
of 3 He along a (100) direction of LiF(OOl). The eikonal 
curves with DFT potential agree with the experimental 
values except for m = at low energies Ei± < 1.0 eV. We 
point out again that the eikonal approximation is a semi- 
classical theory which fails in the vicinity of the classical 
rainbow O r b- The angular position of the quantum rain- 
bow m = is smaller than the classical rainbow angle 
Q r b when X± becomes large 33]. Since the intensity at 
Q r b is overestimated the peak maximum m = is shifted 
to larger deflection angles at lower Ei±. For a correct 
description of the intensity at Q r b a "uniform approx- 
imation" (3ll [32| is necessary, but supernumeraries are 
not affected by this deviation. Since the positions of the 
supernumeraries agree well in both channeling directions 
and since a change of the corrugation of the Hc-LiF(OOl) 
interaction potential by 0.02 A induces a clear shift in 
the position of the supernumeraries [3, [22| , we conclude 
that the DFT He-LiF(OOl) potential is accurate. 

By employing the hard-wall model from Garibaldi et 
al. [33l | valid for sinusoidal small corrugated potential 
surfaces we obtain good agreement as well. Hard wall 
approximations were applied successful for description 
of scattering of He from LiF(OOl) with thermal energies 
0, 0, [H, Due to the large parallel velocity of the 
He projectile the effective potential is averaged along the 
chains of atoms in the beam direction. The potential con- 
tours of the effective potential for scattering along (110) 
for low energies are almost indistinguishable from sine 
functions. For a sinusoidal hard wall, the intensity /„ of 
a Bragg peak of order n is given by: 



FIG. 5: Eikonal results for 3 He atoms scattered from LiF(OOl) 
along the direction (110) for the perpendicular energy Ei± = 
0.5 eV. Dashed (red) line, eikonal differential probability de- 
rived by integrating the initial position over an unit cell; solid 
(blue) line, similar by using an extended integration area, as 
explained in the text. Dotted vertical lines, theoretical peak 
positions based on the Bragg condition (Eq. [7}. 



In 



: ('^[l + cos8„] 



(8) 



with J n being the Bessel function of order n, n = 
arccos-\/l — (nX±/d) 2 the deflection angle of order n, and 
Az the full corrugation of the sinusoidal hard wall, i.e. 
the normal distance between the maximum and the min- 
imum of a equipotential surface. In Ref. it was shown 
that the solution given by Eq. © is in agreement with 
the exact quantum mechanical solution in a wide range 
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FIG. 6: Similar to Fig. [4] but for scattering along the (100) 
direction. 



of Aj_. Since we are interested in the angular positions 
of the supernumerary rainbows we treat n in Eq. J8J as 
£ K and search for the maxima of this oscillating function 
by searching the zeros of the derivation of Eq. ([8]) where 
for Az(Ei±) the effective corrugation of the present DFT 
potential across the (110) channel was used. Results are 
displayed in Fig. 01 showing good agreement with exper- 
iment and with eikonal results. It allows us to conclude 
that the hard wall model is a good approximation here. 

We then apply the hard wall approximation in order 
to describe the different relative intensities of the Bragg 
peaks as shown in Fig. [7] We fit a sum of Lorentzian 
peaks to the diffraction pattern, where the peak positions 
are given by the Bragg relation (Eq. [7]). We assume, 
that the peak width is the same for all diffraction spots 
[3?i |. Resulting relative peak heights are compared with 
the relative intensities given by Eq. ([5]) with A 
fit parameter. The resulting intensity distributions are 
shown as red/gray curves in Fig. [Jj The values for Az 
obtained from the best fit are given in each panel. Since 
Az is almost constant, the differences in the intensities 
are due to the different de Broglie wavelength Aj_ only. 

The effective corrugation Az of the potential across 
the (110) channel deduced from Bragg peak intensities 
for scattering of 3 He and 4 Hc under different angles of 
incidence is plotted as a function of the perpendicular 
energy Ei±_ in Fig. [H] as open symbols. Each data point 
corresponds to a best fit with Eq. (J8]) as in Fig. [Jj Al- 
though the diffraction patterns are different for the He 
isotopes and different angles of incidence, the data are 
well described by Eq. (jHJ) with an almost constant cor- 
rugation Az. We compare the experimentally derived 
values with the effective corrugation of the present DFT 
potential (red/gray full curve in Fig. [5]) defined as the 
normal distance between the maximum and minimum of 



FIG. 7: Angular distributions for intensities projected in- 
side the annulus (cf. Fig.l) for scattering of He atoms from 
LiF(OOl) along <110> under 0, = 0.99 deg. Upper panel: 
Ei = 0.35 keV, middle panel: Ei — 0.50 keV, lower panel: 
Ei = 0.65 keV. Solid curves represent best fits to data by 
sum of peaks with Lorentzian lineshapes. Numbers denote 
diffraction orders. 



FIG. 8: Effective corrugation Az of the potential across the 
(110) channel as a function of the perpendicular energy Ei±. 
Open symbols, experimentally-derived results in the hard 
wall approximation; solid red curve, values derived from the 
present DFT potential; dashed curve, corrugation obtained 
from the interaction potential from Celli et al. ||; dashed 
dotted curve, corrugation from the potential from Celli et al. 
but with the correct He-F~pair potential from Erratum of 
Ahlrichs et al. [H. 
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The small variation of the effective corrugation with 
the perpendicular energy explains why He isotopes with 
the same Aj_ but different perpendicular energies Eij_ 
produce similar diffraction patterns, even though differ- 
ent regions of the surface potential are probed. As ob- 
served in Fig. [9j helium atoms with the same perpendic- 
ular de Broglie wavelength (A^ = 0.13 A) but different 
perpendicular energies (Ei±=1.63 eV and 1.21 eV) show 
identical eikonal distributions, as in the experiment. The 
scattering processes take place indeed at different dis- 
tances to the surface, but since A^ is the same, similar 
interference patterns appear. 



V. CONCLUSIONS 



FIG. 9: Angular distributions for He atoms scattered from 
LiF(OOl) along the direction (110) with \± = 0.13 A. Two 
different isotopes and impact energies are considered. Red 
circles and full curve, experimental data and eikonal results, 
respectively, for 3.5 keV 3 He; blue squares and dashed curve, 
for 2.6 keV 4 He. 



the equipotential surface obtained by averaging the sur- 
face potential along the (110) axial channel. We observe 
that the theoretical curve is close to the experimental 
data. 

Currently, the interaction potential from Celli et al. @ 
is considered as the best available He-LiF(OOl) potential 
[Tol . Hl| . This potential is actually constructed for inter- 
actions at thermal energies but the analytical expression 
can be evaluated also for energies of some eV. The re- 
sulting effective corrugation shows a strong dependence 
on the perpendicular energy Ei± in contrast to the ex- 
perimental values. We conclude that the potential from 
Celli et al. is not adequate for the description of the 
He-LiF interaction in the eV ran ge. This is consistent 
with the results of Ref. [l(| and |Tl| where the exper- 
iment is sensitive to the attractive potential well of the 
planar averaged potential in the meV range only, but not 
to the corrugation of the repulsive part of the interaction 
potential in the eV domain. 

The repulsive part of the Celli potential is based on 
SCF pair potentials for Hc-Li + and Hc-F~ from Ahlrichs 
et al. [33. Due to an error in the calculation of He- 
F~ repulsion [39j |. which is the dominant contribution of 
the repulsive part of the Hc-LiF(OOl) interaction poten- 
tial, the Celli potential has to be corrected. The correct 
He-F~ potential [H, H§] is in good agreement with re- 
cent Hc-F~ ab-initio pair potentials [4(| 5l[ as well as 
the original He-Li + pair potential [38[ agrees with recent 
He-Li + ab-initio potentials Nevertheless, inserting 

the corrected Hc-F - parameters in the expression for the 
Celli potential, the effective corrugation Az(Ei±) devi- 
ates from the experimental results even more. 



We have studied diffraction patterns for swift He atoms 
colliding grazingly with a LiF(OOl) surface in order to 
test the ab-initio surface potential, obtained from DFT 
by making use of the Siesta code. Angular spectra of 
scattered projectiles are obtained with the DFT potential 
from the surface eikonal approximation, which takes into 
account the quantum interference due to the coherent 
superposition of transition amplitudes corresponding to 
different projectile paths with the same deflection angle. 

For incidence along the (110) channel, the dependence 
on the parallel and perpendicular components of the im- 
pact energy was analyzed. It was found that angular 
distributions in terms of the deflection angle 8 are com- 
pletely governed by the perpendicular energy. Diffraction 
spectra as a function of give information about the sur- 
face potential across the incidence channel by means of 
two different mechanisms, supernumerary rainbow and 
Bragg interferences. We have focused here on supernu- 
merary rainbows which are very sensitive to the corru- 
gation of the surface potential within a unit cell. By 
comparison of eikonal angular spectra with experimental 
distributions we concluded that the DFT model provides 
a good description of the surface potential for perpendic- 
ular energies in the range from £^=0,03 eV up to 3 eV. 
The agreement between theoretical and experimental re- 
sults for the intensity near the classical rainbow angle is 
poorer for smaller £/jj_. This deficiency is attributed to 
the range of validity of the semiclassical models, like the 
eikonal approach, which do not include quantum effects 
with respect to the projectile trajectory. 

For scattering along the (110) direction the calculation 
in the hard wall approximation using Eq. j8|) is in good 
agreement with results from simulation with the soft po- 
tential based on the eikonal model. We conclude that in 
this special case the hard wall approximation is good. 

We found that the potential from Celli et al. @] is not 
adequate for the description of the interaction of He with 
LiF(OOl) in the eV energy domain. The DFT potential 
displays an almost constant effective corrugation Az as a 
function of the perpendicular energy, in agreement with 
the experimental values. As a consequence of this nearly 
constant corrugation, He isotopes with the same perpen- 



8 



dicular dc Broglic wavelength but consequently different 
perpendicular energies give rise to the same interference 
patterns. 
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